Heavy-light decay constants with three dynamical flavors 
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We present preliminary results for the heavy-light leptonic decay constants in the presence of three light 
dynamical flavors. We generate dynamical configurations with improved staggered and gauge actions and analyze 
them for heavy- light physics with tadpole improved clover valence quarks. When the scale is set by rUp, we find 
an increase of «23% in fs with three dynamical fiavors over the quenched case. Discretization errors appear to 
be small (^3%) in the quenched case but have not yet been measured in the dynamical case. 



The computation of leptonic decay constants 
for heavy-light mesons plays a key role in the ex- 
traction of CKM matrix elements from experi- 
ment and has been a focus of lattice QCD calcu- 
lations for many years. Here, we report on the 
first attempt to compute such decay constants in 
the presence of Np — 3 light sea quarks {u, d, 
s). For the generation of the dynamical lattices, 
we use the "Asqtad" action Q], which consists of 
order agO^ , improved staggered quarks and or- 
der ala^jO^ Symanzik improved glue 1^,^. This 
action has small discretization errors for many 
quantities |^,^. For the valence quarks, we em- 
ploy a tadpole-improved (Landau link) clover ac- 
tion and the Fermilab formalism . The analysis 
is done in a "partially quenched" manner: the va- 
lence quark masses are extrapolated to their phys- 
ical values with the sea quark masses held fixed. 
The chiral extrapolation of the sea quark masses 
(at fixed lattice spacing) is performed afterwards. 
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Table |l] shows the parameters of our lattices 
and the current state of this project. As we de- 
crease the dynamical quark masses from large 
values in the three-flavor case, we keep the 
three masses degenerate until the physical strange 
quark mass is reached. We then keep the strange 
quark mass fixed as we further decrease m„^d- 

Configurations with varying dynamical quark 
mass (including quenched and two-flavor configu- 
rations) have been matched in lattice spacing us- 
ing the quantity ri which is defined from the 
static quark potential at shorter distance than rg 
10]. From the value ro ~ 0.5 fm we have found 
ri w 0.35 fm §] in full QCD. Note that the er- 
rors in these phenomenological values may be as 
large as 10% 0; we therefore prefer to use mp in 
setting the absolute scale. Details of the lattice 
generation and light quark results are in Ref . (D . 
Our analysis of heavy-light decay constants fol- 
lows Refs. |,|l§; we describe only the important 
differences below. 

On all lattice sets, we compute quark propaga- 
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Table 1 

Status of heavy-light running. We show the nom- 
inal value of a determined from ri ~ 0.35 fm. All 
0.13 fm lattices are 20^ x 64; all 0.09 fm lattices 
are 28^ x 96. "# gen." is the number of configu- 
rations that have been generated; Jb" is the 
number on which fs has been calculated. 
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tors for 5 light and 5 heavy masses. This gives 
good control over both chiral and heavy quark 
interpolations/extrapolations. In Ref. ||^, with 3 
light quark masses, the chiral extrapolation was 
a major source of systematic error. Here, chang- 
ing from linear to quadratic chiral fits of decay 
constants changes the results by <2%. 

The (tadpole improved) perturbative renormal- 
ization of the heavy- light axial current, Z^^^'^, has 
not yet been calculated. At present, we therefore 
proceed as follows at a«0.13 fm: 

(1) Define boosted couphng af from the pla- 
quette following Refs. P,pl|, and then define the 
1-loop coefficient Ca by Z^'^ = 1 + af Ca- 

(2) Fix the scale from nip on each set, with 
valence quarks extrapolated to physical values. 



(3) Fix (a by demanding that / 



quench 
B 



169 



MeV (the MILC continuum-extrapolated result 
|l|). This gives Z^^'^ on aU the 0.13 fm lattices, 
since Ca is independent of the number of dynam- 
ical flavors. 

Note that Ca is dependent on the heavy quark 
mass through the dimensionless quantity am. To 
test scaling we first repeat the above procedure 



for /c at 0.13 fm and then interpolate to the cor- 
rect quark mass for fs at 0.09 fm. There is a 2% 
change in Ca for fs between the two lattice spac- 
ings. Because this change results in only a 0.3% 
effect on /b, we believe that systematic effects in 
the procedure are negligible. 

Note also that comparisons of Jb with 
or changes of /b with a are only very weakly de- 
pendent on our assumption about the continuum 
value of Further, the ratio /bJIb is es- 

sentially independent of perturbation theory. 

In the plots below, we show two sets of errors. 
The smaller ones represent the statistical errors 
from the fits, computed by jackknife. We use the 
same fitting ranges in time for all lattices with 
the same spacing. The larger error bars take into 
account variations caused by changing the fitting 
ranges while keeping confidence levels and statis- 
tical errors reasonable. 

Preliminary data for /b with the scale set by 
nip is shown in Fig. |l} There is clear evidence for 
an increase in /b as the dynamical quark mass is 
decreased. There is also an increase in the Np = 3 
results over those with Np = 2. Note that the 
quenched results at the two lattice spacings are 
consistent; they differ by <3%. 
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Figure 1. /b versus (m^/mp)^ of the dynamical 
light quarks. The scale is set by mp, with valence 
quarks extrapolated to the physical point. 

Figure || shows /b with the scale set by ri = 
0.35 fm (the phenomenological error in this value 
of ri is not included). There is now little if any 
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dependence on the dynamical quark mass: Hold- 
ing ri fixed forces the short distance potentials to 
match quite closely, which in turn keeps /b, the 
"wave function at the origin," more or less fixed. 
In the physical case (three dynamical quarks ex- 
trapolated to physical masses) the values of fs in 
Figs. |l| and I are consistent. 
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Figure 2. Same as Fig. |l|, but with scale set by 
ri = 0.35 fm. 

The ratio fs^ / fs is shown in Fig. ||. No obvi- 
ous dependence on the number or masses of dy- 
namical quarks is apparent. We therefore fit all 
the data to a constant. 
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Figure 3. Same as Fig. [|, but for /bJ/b- 

For the moment, we assume that the scaling 
errors in Jb or /b,/ fs arc no larger than the dif- 



ference in our quenched results at 0.13 fm and 

0. 09 fm. With this assumption (which will be 
tested in the coming year) we arrive at the follow- 
ing preL'minary results: — 1.23(4)(6) 
and fB,/ fs = l-18(l)(ii), where the first error 
is statistical and the second is systematic, includ- 
ing discretization and chiral extrapolation (va- 
lence and dynamical) errors, the uncertainties in 
Kg (for JbJ Is), and a rough estimate of pertur- 
bative errors. A direct determination of Jb itself 
with three dynamical flavors must await the cal- 
culation of the renormalization factors. 
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where /,r sets the scale for the central value. 
With an mp scale, this becomes 169 MeV. 



